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In this work we study the contribution to leptogenesis from AL = 1 decay and scattering pro- 
cesses mediated by the Higgs with quarks in the initial and final states using the formalism of 
non-equilibrium quantum field theory. Starting from fundamental equations for correlators of the 
quantum fields we derive quantum-corrected Boltzmann and rate equations for the total lepton 
asymmetry improved in that they include quantum-statistical effects and medium corrections to 
the quasiparticle properties. To compute the collision term we take into account one- and two-loop 
contributions to the lepton self-energy and use the extended quasiparticle approximation for the 
Higgs two-point function. The resulting CP-violating and washout reaction densities are numeri- 
cally compared to the conventional ones. 
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I. INTRODUCTION 

The Standard Model (SM) of particle physics [iHll 
has successfully passed the numerous experimental tests 
performed so far. The recent observation of the Higgs 
particle Q at the LHC [a, Q also seems to confirm the 
mechanism of spontaneous symmetry breaking, which is 
responsible for masses of the known gauge bosons and 
fcrmions. On the other hand, wc know that the SM is 
not complete. Firstly, it does not provide a viable dark 
matter candidate. Secondly, it predicts that the active 
neutrinos are strictly massless, which contradicts the re- 
sults of neutrino oscillation experiments. A simple yet el- 
egant way to generate small but nonzero neutrino masses 
is to add three right-handed Majorana neutrinos to the 
model: 

- h,JJPRN, - hi^N,4i^ PlIo. , (1) 

where Nt = iVf are the heavy Majorana fields, ia are the 
lepton doublets and (j) = ia2(t>* is the conjugate of the 
Higgs doublet. After the electroweak symmetry break- 
ing the active neutrinos receive naturally small masses 
through the type-I seesaw mechanism. This scenario has 
even more far-reaching consequences as it can explain 
another beyond-the-SM observation, the baryon asym- 
metry of the universe. The Majorana mass term in ^ 
violates lepton number. In the early Universe a decay 
of the Majorana neutrino into a Icpton-Higgs pair in- 
creases the total lepton number of the Universe by one 
unit, and a decay into the corresponding antiparticles 
decreases the total lepton number by one unit. If there 



is CP-violation then, on average, the number of leptons 
produced in those decays is not equal to the number of 
antileptons and a net lepton asymmetry is produced. It 
is also known that whereas the difference of the lepton 
and baryon numbers is conserved in the Standard Model, 
any other their linear combination is not [7[ . This implies 
that the lepton asymmetry produced by the Majorana 
neutrinos is partially converted to the baryon asymmetry 
[g . This mechanism, which is referred to as baryogenesis 
via leptogenesis, naturally explains the observed baryon 
asymmetry of the Universe. For a more detailed review 
of leptogenesis see e.g. j9l-[lll|. 

The state-of-the-art analysis of the asymmetry gener- 
ation uses Boltzmann equations with the decay and scat- 
tering amplitudes calculated in vacuum. Their applica- 
bility in the hot and expanding early universe is ques- 
tionable and can be cross-checked using a first-principle 
approach based on the use of non-equilibrium quantum 
field theory. One of the most important processes for the 
generation of the asymmetry is the decay of the Majorana 
neutrino. Thermal effects enhancing CP- violation in the 
decay have been studied in [12 -llal . The role of the flavor 
effects has been addressed in 13]. A first-principle anal- 
ysis of the asymmetry generation in the very interes ting 
regime of resonant leptogenesis has been presented in [l8| 
and [13 ■ The effect of next-to- leading order corrections 
from the gauge interactions of lepton and Higgs doublets 
on the production and decay rate of right-handed neu- 
trinos at finite temperature has been recently studied in 
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The asymmetry generated in the Majorana decay is 
partially washed out by the inverse decay and scattering 
processes. The latter can be classified into two categories. 
The first category includes AL ~ 2 scattering processes 
mediated by the Majorana neutrinos. A first-principle 
analysis of such processes free of the notorious double- 
counting problem has been presented in j22| . The second 
category includes AL = 1 decay and scattering processes 
mediated by the Higgs. The latter processes are also 



known to play an important role in the asymmetry gen- 
eration and are addressed in the present paper. 

The outline of the paper is as follows. In Sec.|ll] we 
briefly review the canonical approach to the analysis of 
the AL = 1 processes and derive the corresponding am- 
plitudes and reduced cross-sections. In Sec. lIIII we derive 
quantum-generalized Boltzmann equations for the lep- 
ton asymmetry, calculate the effective amplitudes of the 
Higgs-mediated scattering processes and compare them 
with the canonical ones. The obtained Boltzmann equa- 
tions are used in Sec.|V] to derive a simple system of 
rate equations for the total lepton asymmetry. In sec- 
tion Scc. lVIl we present a numerical comparison of the 
corresponding reaction densities with the ones obtained 
using the canonical approach. A summary of the results 



is presented in Sec. lVII 



II. CONVENTIONAL APPROACH 

In the scenario of thermal leptogenesis lepton asymme- 
try is generated in the lepton number and CF-violating 
decay of the heavy Majorana neutrinos. The corrcspond- 
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FIG. 1: Tree-level, one-loop self-energy and one-loop vertex 
contributions to the decay of the heavy Majorana neutrino. 

ing CP-violating parameters receive contributions from 
the interference of the tree-level amplitude with the ver- 
tex 0, [l^ and self-energy [ij, |23| - [26| amplitudes, see 
Fig-Hi The contribution of the loop diagrams can be ac- 
counted for by effective Yukawa couplings [2g . If thermal 
masses of the SM particles are negligible, they are given 
by: 



"'+,0:1 — f^ai '^i^ajyi^ '^Jji Qij 

where the loop-function gij is defined as 
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Note that this expression is valid only for on-shell final 
states. The first term in ([3]) is related to the self-energy 
and the second term to the vertex contribution. This 




FIG. 2: Effective one-loop diagram for the self-energy and 
vertex contributions to the decay of the lightest Majorana 
neutrino for a strongly hierarchical mass spectrum. 



expression is applicable for a mildly or strongly hier- 
archical mass spectrum of the Majorana neutrinos. In 
both cases most of the asymmetry is typically generated 
by the lightest Majorana neutrino, whereas the asymme- 
try generated by the heavier ones is almost completely 
washed out. For a strongly hierarchical mass spectrum, 
Mi <C Mj, the intermediate Majorana line in Figs.[T]b 
and[T]c contracts to a point, see Fig. [51 and the structure 
of the self-energy and vertex contributions is the same. 
In this limit: 



3_m 



(4) 



Note that in this approximation the loop integral leading 
to ^ depends only on the momentum of the initial state 
and is independent of the momenta of the final states. 
This implies in particular that this expression can also 
be used for off-shell final states. 

Using the effective couplings ^ we find for the decay 
amplitudes (squared) |22l.[26|: 

SjVi^f0 = 9N9w{h\h^)i,{qp) , (5a) 

-Ni^i4, = 9N9w{h''-h_)ii{qp) , (5b) 

where we have summed over fiavors of the Icptons in the 
final state as well as over the Majorana spin (g^ ~ 2) and 
the SU{2)l group (g^, = 2) degrees of freedom. Here q 
and p are momenta of the heavy neutrino and lepton, 
respectively. The decay amplitudes (O can be traded for 
the total decay amplitude and CP-violating parameter: 






(6a) 
(6b) 



Combining (O and ^ we then find for the (unflavored) 
CF-violating parameter: 



Im(fetfe)2^. 



X 25ij , j ^i. 



(7) 



The asymmetry generated by the Majorana decay is par- 
tially washed out by the inverse decay and scattering 
processes violating lepton number. An important role is 
played by the AL = 2 scattering processes mediated by 
the heavy Majorana neutrinos [22, lla, [27| . In addition, 
there are AL = 1 scattering process mediated by the 



Higgs doublet with quarks (or the gauge bosons) in the 
initial and final states [Mllj, see Fig.Eland Fig.H The 
Higgs coupling to the top is considerably larger than to 
the other quarks of the three generations. For this rea- 
son we do not consider the latter here. The corresponding 
Lagrangian reads: 



■^SM D ~\Qc^PRt - \*iPL(f>^Q , 



(8) 



where Q and t arc the SU{2)l doublet and singlet of 
the third quark generation. The AL = 1 processes are 




FIG. 3: Tree-level, self-energy and vertex contributions to 
the scattering processes NiQ — ^ It. Similar diagrams for the 
scattering process Nit —J- £Q are obtained by replacing Q with 
i and t with Q as well as inverting the direction of the arrows. 

also CP-violating. The CP-violation is generated by the 
same self-energy and vertex diagrams. Strictly speaking, 
since the Higgs is no longer on-shell the effective cou- 
plings ^ arc not applicable in this case. On the other 
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where Aj,(fc) 



i/{e 



is the Feynman (or time- 



ordered) propagator^ of the intermediate Higgs and we 
have defined 



"0Q^t =25s|A|^(pQPt)- 



(10) 



Here gs = 3 is the SU{2i)c factor, and pt and pq are 
the momenta of the singlet and the doublet respectively. 
For the charge-conjugate process we find an expression 
similar to ©. As can be inferred from (fTU]) in this work 
we neglect CP-violation in the quark sector, which is 
known to be small. Defining CP-violating parameter in 
scattering as 



ejf^Y = — —^ 



(11) 



we then obtain for eNiQ^it the same expression as for 
the Majorana decay, see ([7]). In the same approximation 
amplitude and CP-violating parameter for Nit — > iQ 
scattering coincide with those for NiQ — > U process. 
Proceeding in a similar way we find for the scattering 
amplitude of the Nil —>■ Qt process: 



^Nil^Qt ~~ ^Nit- 



xAl{p + q) XS. 



'rQt 



(12) 



where S .^^^ = :zj^q^^ because we neglect CP-violation 
in the quark sector. Furthermore, for a strongly hierar- 
chical mass spectrum S^^ , ~ Sjy._^^i,. The resulting 
expression for the CP-violating parameter then coincides 
with ©. 

If the lepton and both quarks are in the final state 
then instead of a scattering process we deal with a three- 
body Majorana decay, see Fig.[Sl In complete analogy 
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FIG. 4: Tree-level, self-energy and vertex contributions to the 
scattering processes Nil -^ Qt. 

hand, for a strongly hierarchical mass spectrum the in- 
termediate Majorana lines in Fig. [3] and Fig. [5] again con- 
tract to a point and the momenta of the Higgs and lepton 
play no role. In other words, for a strongly hierarchical 
mass spectrum we still can use the effective couplings 
([2|) supplemented with Q to calculate the CP-violating 
scattering amplitudes. 

Summing over flavors and colors of the quarks and lep- 
tons in the initial and final states as well as over the cor- 
responding SU{2)l and spin degrees of freedom we find 
for the amplitude of NiQ — > £t scattering: 
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(9) 
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FIG. 5: Tree-level, self-energy and vertex contributions to the 
amplitude of the three-body decay processes Ni -^ IQt. 

with the scattering processes we can write its amplitude 
in the form: 



^-N.^iQt ~ ^N.^e<p ^ ^t{pq +Pt) X ^4,^Qt ■ 



(13) 



Evidently, CP-violating parameter for this process coin- 
cides with that for the two-body Majorana decay. 



^ In the kinematic region of interest the decay width term in the 
Feynman propagator of the Higgs plays no role and can bo ne- 
glected. 



To compute the generated Icpton asymmetry, the con- 
ventional approach uses the generalized Boltzmann equa- 
tion for the total lepton abundance, Yt = nj^j s. with 
s being the comoving entropy density [28|. In the 
Friedmann-Robertson- Walker (FRW) universe the con- 
tribution of the Higgs-mediated processes to the right- 
hand side of the the Boltzmann equation simplifies to: 

sUdYL _ 

z dz 

-E 



E 



Na^Qifrnfii'^ - /q)(1 - /t) - inverse] 

[Sjv.q^^Jat./qCI - /<?)(! - ft) - inverse] 
a''^%l^i2nrS{q + pt-p-pQ) 

^ [-Nj^iQfNjiC^ -/£)(!- /cj) - inverse] 

+ E / '^KaQt' (2-)''5(<z -P-PQ-Pt) 
i •' 

^ ["w.->«Qt/jv.(l - /£)(! - /q)(1 - ft) ~ inverse] 

— CP conjugate processes. (14) 

where we have introduced the dimensionless inverse tem- 
perature z = Mi/T, the Hubble rate % = H\t=Mi, 
and dn^^^'j,'^"''^ " stands for the product of the invari- 
ant phase space elements, dng = d^p/\{2'KY'2Ep\. Note 
that to ensure vanishing of the asymmetry in thermal 
equilibrium one should also include CP-violating 2 o 3 
processes [l^l- Since there is no need for that in the 
non-equilibrium QFT approach we will not consider these 
processes here. 



III. NON-EQUILIBRIUM QFT APPROACH 

The formalism of non-equilibrium quantum field the- 
ory provides a powerful tool for the description of out-of- 
equilibrium quantum fields and is therefore well suited for 
the analysis of leptogenesis. In this section we briefiy re- 
view results obtained in [22l and introduce notation that 
will be used in the rest of the paper. As has been argued 
in [221, the equation of motion for the lepton asymme- 
try can be derived by considering the divergence of the 
lepton current. In the FRW Universe j^ = (nL,0) and 
therefore it is related to the total lepton abundance by: 






sHdYL 

z dz 



(15) 



Using the formalism of non-equilibrium quantum field 
theory one can express it through propagators and self- 
energies of leptons. After a transformation to the Wigner 



space we obtain [22l |: 

V^j^it.p) ^ g^ J dnltT[t^it,p)S^it,p) 

-S>(t,p)5<(t,p)], (16) 

where dUi = (i^p/(27r)* and the hats denote matrices 
in flavor space. In p6p we have taken into account that 
the SU{2)l symmetry is unbroken at the epoch of lep- 
togenesis. As a consequence, the SU{2)l structure of 
the propagator is trivial, 5""^ = SabS"^ , and summa- 
tion over the SU{2)l components simply results in the 
overall factor g^ = 2. Furthermore, in this work we re- 
strict ourselves to the analysis of unfiavored leptogenesis. 
Therefore, the lepton propagator can be approximated by 
S'"'^ = 5"^ S . Similar relation also holds for the lepton 
self-energy. Then the equation for the divergence of the 
lepton current takes the form: 



v„i>i 



I^Jl — 9n 



dp^ f d^p 



27r J (27r)3 



(17) 



tr[(S<%-E>^<)-(S<5>~E>5<)] 



where E = S"" and we have suppressed the argument 
{t,p) of the two-point functions. Note that the trace 
in (|17p acts now in spinor space only. To convert the 
integration over positive and negative frequencies into 
the integration over positive frequencies only we have in- 
troduced in (J17p CP-conjugate two-point functions and 
self-energies which are denoted by the bar. Accordin g to 
the extended quasiparticle approximation (eQP) J29l - l3l| 
the Wigthmann propagators can be split into on- and 
off-shell parts: 



S'> — S'> — 2('S'i^'S>5^ + Sj^Y.^Sji 



(18) 



The off-shell parts of the lepton propagators exactly can- 
cel out in the lepton current as they are lepton number 
conserving. On the other hand, as we will see later, the 
off-shell part of the Higgs two-point functions is crucial 
for a correct description of the scattering processes. The 
on-shell part of the Wightman propagators is related to 
the eQP spectral function and one-particle distribution 
function /^ by the Kadanoff-Baym (KB) ansatz: 



S*. 



(1 - fi) Sp , S'< 



-ft Sp , 



where 



Sp — ^i^^r'^p^r'^pSa^pSa ■ 



(19) 



(20) 



In the limit of vanishing width the eQP spectral function 
Sp approaches the Dirac delta- function |22| . 

Sp « (2^) sign(/)5(p2 - ml)PLfPR 

^Sp{p)PLfPR, (21) 



where we have extracted the 'scalar' part S for nota- 
tional convenience. In (|2ip we have approximately taken 
the gauge interactions into account in the form of effec- 
tive masses of the leptons. Note that we will not attempt 
a fully consistent inclusion of the gauge interactions here. 
In the used approximation the spectral function is CP- 
symmetric. This implies that the spectral properties, in 
particular the masses, of the particles and antiparticles 
are the same. 

To evaluate the right-hand side of (fT7)) we need to spec- 
ify the form of the lepton self-energy. It can be obtained 
by functional differentiation of the 2PI effective action 
with respect to the lepton propagator. Loosely speaking, 
this means that the self-energies are obtained by cutting 
one line of the 2PI contributions to the effective action. 
The two- and three-loop contributions are presented in 
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FIG. 6: Two- and three-loop contributions to the 2PI effective 
action and the corresponding contributions to the lepton self- 
energy. 

Fig-ISa) and Fig.EI^c). The one- loop contribution to the 
lepton self-energy, see Fig. [6] (b) . is given by [22| : 

^%\t,p) = - j <m%{2n)H{p + k-q) 

X {h^h),,PR.y^{t,q)PL/^^{t,k), (22) 

where 5^ and A denote the Majorana and Higgs propa- 
gators respectively, and dll^j. = dll^ dll^. . The expres- 
sion for the two-loop contribution, sec Fig.[S](d), is rather 
lengthy. Here we will only need a part of it: 

Y.f {t, p) = J dn^, {27T)^S{p + k-q) (23) 

X [{h^hh^{h^h)j,nA^^{t,q,k)PLCS^^{t,q)PLA^{t,k) 
+ ih^hU{h^h)„,,PR.y^'{t,q)CPRVnrn{t,q,k)A^{t,k)] , 

where we have introduced two functions containing loop 



Aran{t,q.k) = j dlli^ 



fc2fc3 



X {2Ti)H{q + fci + fca) {2TTf5{k + k2~ fcs) 
X [PR^rit, -k^)CPRSl{t, k2)AA{t, fci) 
+ PR^rii. -k3)CPRSl{t, fc2)A^(t, fci) 
+ PRS^rit,-k3)CPRS^{t,k2)Ap{t,ki)] , 



(24) 



and Vnm{t,q,k) = P Al^„^{t,q,k) P to shorten the nota- 
tion. Here P = 7^^ is the parity conjugation operator. 
The remaining terms of the two-loop self-energy can be 
found in |22| . As has been demonstrated in the same ref- 
erence, CP-conjugates of the above self-energies can be 
obtained by replacing the Yukawa couplings by the com- 
plex conjugated ones and the propagators by the CP- 
conjugated ones. 

Comparing ([^^ and ([^5]) we see that the two self- 
energies have a very similar structure. First, the integra- 
tion is over momenta of the Higgs and Majorana neutrino 
and the delta-function contains the same combination of 
the momenta. Second, they both include one Wightman 
propagator of the Higgs field and one Wightman prop- 
agator of the Majorana field. These can be interpreted 
as cut-propagators which describe on-shell particles cre- 
ated from or absorbed by the plasma [321 • The retarded 
and advanced propagators can be associated with the 
off-shell intermediate states. We therefore conclude that 
the two self-energies describe CP-violating decay of the 
heavy neutrino into a lepton-Higgs pair. Note that this 
interpretation only holds for the "particle" part of the 
eQP ansatz. The inclusion of the off-shell part of the 
Higgs Wightman propagator gives raise to the Higgs me- 
diated scattering processes and three-body decay, see sec- 
tion |TVl 

To evaluate ([^^ and (P5|) we need to know the form of 
the Higgs and Majorana propagators. For the Higgs field 
we will adopt in this section a leading-order approxima- 
tion: 



A> = (1 + /^)A^ 



A, 



U\: 



(25) 



and a simple quasiparticlc approximation for the spectral 
function, 



A^(i,fc)= (27r)sign(fcO)^(fc2-m^)^ 



(26) 



where m^ is the effective thermal Higgs mass. Close to 
thermal equilibrium the full resummed Majorana propa- 
gator is given by [23 |: 



^^ 



-^&r[S^ 
SrU^S^ - i(4ll|4 + 5^Il|5^)]e^, (27) 



where H'' and H denote the diagonal and off-diagonal 
components of the Majorana self-energy respectively, 5^ 
and 5^ are given by 

-1 



S 



R{A) 



M 



^R{A)] 



(28) 



and wc have introduced 
-1 



^H^ 



1- 



SrH^) 



Qa^ 



1- 



^aSa) 



(29) 



to shorten the notation. The first term in the square 
brackets of l^7}i describes (inverse) decay of the Majo- 
rana neutrino, whereas the remaining three terms de- 
scribe two-body scattering processes mediated by the 
Majorana neutrino. For tlie "particle" part of the eQP 
diagonal Wightman propagators of the Majorana neu- 
trino one can use the KB approximation: 



C7in 



{i-fNjs;' 



onn 



~fN,.s: 



(30) 



with the spectral function given by an expression identi- 
cal to (j20p . Substituting p8)) we find in the limit of small 
decay width: 






(2^)sign(g")%2_M2)(^ + M„) 



^S^ (^ + M„). (31) 

Inserting (P^ and ([221) into the divergence of the lepton 
current (|17p and integrating over the frequencies we then 
obtain an expression that strongly resembles the Boltz- 
mann equation: 



sHdYL 

z dz 



-E 



JTT9Pfe 



X 






q-pk 



where we have introduced 



__-pq;pk _] 



(32) 



-r-PaPb--lPiPj-- _ /'0^^4e/ 
^ab..^ij.. = (27r) d{Pa 



[f^^f!' 



lb 
fPj 



Pb + ■ ■ ■ - Pi - Pj - 



IT) 



/f/f ...(i±/r)(i±/r)---] 



(33) 



with the plus (minus) sign corresponding to bosons 
(fermions). Note that ^^^^^j'^'^^" vanishes in equilib- 
rium due to detailed balance. This implies that in accor- 
dance with the third Sakharov condition [33| no asymme- 
try is generated in equilibrium. In the Kadanoff-Baym 
formalism this result is obtained automatically and no 
need for the real intermediate state subtraction arises. 

The effective decay amplitudes S are given by a sum 
of the tree- level, one- loop self-energy and one-loop vertex 
contributions. The first two: 












X tv[ei^ 



9w2^mn\'^ '''Iran 

X tr[e^'(g)(^ + M{)QT{<l)PLfPE] 



(34a) 



7:3 



(34b) 

emerge from the one- loop lepton self-energy (1^^ . The 
third one: 






t:V 



^^ = ~ gn,{h^h)l M, tr[%(<z, k)CPLfPR] 
- gUh^h)% M,tT[CVjjiq, k)PLfPR\ , 

~ gUh^h)% M, tr[%((Z, k)CPL^PR\ , 



(35b) 



is generated by the two-loop lepton self-energy ([^B]) . Sub- 
stituting ([M)) and (|551) into (O we find to leading order 
in the couplings that the total decay amplitude summed 
over the Majorana spin degrees of freedom is given by 
S^. = 2qNQw {h^h)ii (pq) ■ The self-energy CF- violating 
parameter reads J22l | : 



E 



lmih^h)l NhVjpLs 



qp 



{h^h)u{hm)jj M] 



M^S'^iq), (36) 



where the 'scalar' part of the diagonal hermitian Majo- 
rana propagator is given by [22| : 



'5i^('z) 



h[S'^il) 



q^ - i\/| 



{q' 



M2)2 



{r,/M, ■ qLs 



(37) 



It describes the intermediate Majorana neutrino line in 
Fig.[T]b. Note that ((5^ has been obtained assuming a 
hierarchical mass spectrum of the heavy neutrinos and 
is not applicable for a quasidegenerate spectrum. For 
positive g^ and q^ the self-energy loop function Ls is 
given by j22| : 

L^ = WttJ dUt1.^^i27rr5{q ~k,~p,)p'^ 

x[i+4^-/r]- (38) 

Simplifying psp we find for the vertex CF- violating pa- 
rameter [22|: 

v_ ly Im(fe^fe)fj M,T,pLv 

'^ 2 2^Q,muhth),, M| qp ■ ^^""^ 

The vertex loop function is given by: 

L^{q,p)^16TTMfl dn*^^^^^ (40) 

X {2n)Hiq + ki+pi){2Tr)^S{q-p + pi-qi)p'^ 



[A^{k,)SAPi)Si'{qi) 

A,(fci)S,(pi)5].^(9i) + 

.A,(fci)S,(pi)5].^(gi) + 



FA^(fci)S,(pi)5^,^(gi) 

A^ik,)SAPi)S'p'{qi) 
AAki)S,,{p,)S^^^q,)] , 



where Sp = (<S^ + S^)/2 is the 'scalar' part of the cor- 
responding statistical propagator of the heavy neutrino. 
For the lepton and Higgs fields the definitions are similar. 
The three lines in the square brackets in p9p correspond 
to different cuts through two of the three internal lines of 
the vertex loop. The first line corresponds to cutting the 
propagators of the Higgs and lepton and can be simplified 
to M: 



(35a) p4^9.P) = 16^/ dntl^i2Tr)^S{q-pi - k,) 



x(m)[i + 4^-/r] 



M^ 



Mf-iq-pi-pJ^ 



(41) 



The other two are cuts through the Majorana and lepton 
and the Majorana and Higgs Uncs respectively [15|- For 
the second cut we obtain: 



-Nil, 



pLy^'iq,p) = 16TT / rfn^%^^ {2ny6iq - p + Pi - qi) 



Mi 



' TO^- {q + piY 

+ 167r I <m^^ dlil^ {2^f5{q -p-pi+qi) 



x(m)[/l?.-/r] 



A/| 



rnl- [q-piY ' 
whereas contribution of the third cut is given by: 



(42) 



-N,<p, 



N,4, 



pLy'''{q,p) = 16^ / dn;;^(2^)^J(qi -p- fci) 



x(p<Z+pfci)[4^+/^.] 



M] 



m'j - {q + kiY 



, (43) 



where we have assumed Mi < Mj so that the (inverse) 
decay Ni o Nji£ is kinematically forbidden. In P^ the 
second term vanishes for the decay process Ni o £(j) but 
gives a non-zero contribution for the scattering processes, 
see section ITVl If the intermediate Majorana neutrino is 
much heavier than the decaying one the last two cuts are 
strongly Boltzmann-suppressed. Furthermore, compar- 
ing ([55| and ([iT|) we observe that in this case pLy ~ pLs- 
In the same approximation we can also neglect the 'reg- 
ulator' term in the denominator of p7|) . The two con- 
tributions to the CP-violating parameter then have the 
same structure and their sum can be written in the form: 



qp 



(44) 



In the vacuum limit Lg = q^ and we recover ([7|) . At finite 
temperatures the CP-violating parameter is moderately 
enhanced by the medium effects [22 1. 



IV. HIGGS MEDIATED SCATTERING 

In the previous section we have approximated the full 
rcsummcd Higgs propagator by leading-order expressions 
(PSj) and (|26p . In this section we will use a more accurate 
eQP approximation. As we will see, it allows one to de- 
scribe Higgs-mediated AL = 1 two-body scattering and 
three-body decay processes. 

Similarly to (jlSp . the extended quasiparticle approxi- 
mation for the Higgs propagator reads: 



^^-h-H^R 



Al)fi^ 



(45) 



Its graphic interpretation is presented in Fig.[71 For the 
first term on the right-hand side of ((45|) we can again use 
approximations ([^S]) and (|26|) . To analyze the second 



(a) 



(b) 



(c) 



(d) 



FIG. 7; Schematic representation of the eQP approximation 
for the Higgs field. 



term we have to specify the Higgs self-energy. At one- 
loop level it reads: 



n^ (i, k) =g,\X\' J dH^^p^ (2^)4<5(fc -pt+ pq) 

X tr[SQ^{t,pQ)PRSt;,it,pt)PL] , (46) 

see AppcndixEl for more details. As is evident from (|46p . 
here we limit our analysis to contributions generated by 
the quarks of the third generations. Let us note that in 
the SM the gauge contribution to the Higgs self-energy 
is of the same order of magnitude and should not be 
dismissed in a fully consistent approximation. Using the 
KB ansatz for the cQP propagators of the quarks with 
effective thermal mass: 



St>^il-ft)Stp, St^ = -ftStp, (47a) 

^Q> = (1 ^ fQ)SQp , Sq^ == ^IqSqp , (47b) 

~Sqp = (27r)sign(pc3°)<S(PQ' - uiDPlP^Pr^ 

= SQ^PLfQPR, (48a) 

Stp - {27r)sign{pt°)SipY - mf)PRf^PL 

= StpPRf^PL , (48b) 



with 



and neglecting their off-shell parts, which are lepton num- 
ber conserving, we can write the Higgs self-energy in the 
form: 

r!> (t, k) = -2.9, 1 Ap ydH^^p, (27r)4<5(fc + pQ - p^) 

X /q(1 - ft){PQPt)SQpipQ)Stp{pt) , (49a) 

n^{t,k) = -2g,\X\'JdU^^^j,^{27r)H{k + pQ-pt) 

X (1 - fQ)ftiPQPt)SQp{pQ)Stp{pt) , (49b) 

Substituting the one-loop lepton self-energy ((22|) with the 
Higgs propagator given by (gS]) into the divergence of the 
lepton current (|17p . we obtain: 



sHdYL 

z dz 



= ^ / dH^pQpp, (27r)^(5((7 +PQ-p-pt) 



xS^iq)S^{p)SQ^{pQ)Stpipt) 

X Ejy^^g^{q,p)AJi^^{pt - PQ)E^Q^t{pt,PQ) 

x[/^C(i-/r)(i-/r) 

-/r/r(i-/^.)(i-C)], (50) 
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where we have introduced a combination of the retarded 
and advanced propagators, 



K+Aik) 



i[A?,Xfc)+A^(fc)] 



(51) 



which describes the intermediate Higgs hne in Fig. [2] and 
Fig.m Note that in ((50)) the momenta are not restricted 
to the mass shell. In particular, the zeroth components 
of the momenta can have either sign. Due to the Dirac- 
deltas in the spectral functions the frequency integra- 
tion is trivial. Each of the spectral functions can be de- 
composed into a sum of two delta-functions, one with 
positive and one with negative frequency, leading to 2^ 
terms. These different terms correspond to 1 -H- 3 (in- 
verse) decay, 2 o 2 scattering and to (unphysical) -(-> 4 
process. An additional constraint comes from the delta- 
function ensuring energy conservation. In the regime 
Mi > mi + mq + mi only 8 terms satisfy the energy 
conservation. Using the relation 



l±/„(t,-p) = T/a(i,P): 



(52) 



where fa denotes the distribution function of the antipar- 
ticles, we can then recast ([50)1 in the form: 



sUcLYl 

z dz 



■E/^n-7, 



pt 

iQt 



(53) 



-pIPQ 'PPt ■=■ 



-pQPQJPPt ^ _ _ 



-pqpt\PPQ_ 77 -pqpt \PPQ .^ 



■pIP'PQPtT: 
-pQ-PPQPt .^ 



-pQp;PQP_t ■^ _ _ 

-pl^PPQPt ^ 

•^ Ni-^lQt^Ni^£Qt 



The effective scattering amplitudes in (|53)l correspond to 
different assignments for the sign of the four-momenta 
in (|50p . reflecting the usual crossing symmetry. For the 
tree-level and self-energy contributions to the effective 
scattering and decay amplitudes we obtain: 






-T+_S 
■:T+S 



^T+S a2 / 

^T+S a2 / 

^T+S a2 / 
■^N,^e4>^R+A\Pt 



' PQ)^4>^Qt ' 



(54a) 
(54b) 
(54c) 



and similar expressions for the CP-conjugate ones. Note 
that SjY^^^ J and S^^^^ are given by the same expres- 
sion since the CP-violating loop term in (|34[1 depends 
only on the momentum of the Majorana neutrino. In 
vacuum these scattering amplitudes reduce to (O and 
(J12p respectively but with the Feynman propagator A|n 
replaced by A^_,_^ . In the latter the contribution of the 
real intermediate state is subtracted by construction J22l | . 
However, in the regime m^ < m,Q + m,t the intermediate 
Higgs cannot be on-shell such that the vacuum and in- 
medium amplitudes become numerically equal. Since the 



amplitudes S ,^^^ and S ^q^^ factorize in (|5^ and are 
CP-conserving, the self-energy CP-violating parameter 
in these processes is the same as in the Majorana decay, 
see ([M]) . However, since the decay and scattering pro- 
cesses have different kinematics the averaged decay and 
scattering CP-violating parameters are not identical. 

Next we consider the two- loop lepton self-energy (|23p . 
Proceeding as above we find for the divergence of the 
lepton current an expression of the form (j53p with the 
amplitudes given by: 






~V _ _ 






[Pt 



^R+A\Pt 

'-'Nd^ij,^R+A{Pt 
'-'Ni'^irt>^R.+A\Pt 



PQ)^4>Q<^t J 

PQ)^(f>^Qt ' 

' PQJ^dx^Qt ■ 



(55a) 
(55b) 
(55c) 



Since the vertex contribution to the Majorana decay am- 
plitude depends on the momentum of the Higgs, the am- 
plitude 5](^.^^i, does not coincide with S]^^^ , and we 
can define two inequivalent vertex CP-violating param- 
eters [3J|- For the scattering processes NiQ o it and 
Nit O £Q as well as for the three-body decay Ni -(->• iQt, 
the CP-violating parameter coincides with p9p with the 
contributions of the three possible cuts given by (|4ip . 
dH]) and (03]) respectively. For the N^i <-^ Qt process, 
the CP-violating parameter still has the form ((39)) . but 
since the lepton is in the initial state the loop integral 
must be evaluated at {q,—p) instead of {q,p). For the 
first cut we obtain: 

pL'y^{q,p) = 16^1 dHf^p^(2^)4^(g-pi - fci) 



fPil 



Mf 



xta,)[l + /?-/f,,„_,,_^,^^,. 



(56) 



Contributions of the second and third cuts are given by: 
pLp'{q,p) = IGttJ d\il^^^{2TTf5{q+p-pi - qi) 

X (pp.) [1 - e - /r] ^ J 

- 16^ / dn^^^l {2T,)H{q +p + pi-qi) 



x(ppi)[/^^-/r](^ 



Mf 



+ pi)2 -ml 



(57) 



and by 

pL^^^q^p) = WttJ dH^'^(2^)45(qi -p-h) 

x(M^.fco[4-+/^.] (^,,^;.,^. , (58) 



respectively. As follows from ([55]) and ([57)) . CP-violating 
parameter in the N ji -o - Qt scattering receives two vac- 
uum contributions [34j . One is the usual cut through 



l(j), and the second one is given by the first term in the 
cut through Nji. The kinematics of the second contri- 
bution corresponds to Nil <-!• Njl t-channel scattering 
and therefore requires the initial center-of-mass energy 
s = g + p to be greater than the final masses Mj + tti^, 
meaning that contribution of this term to the reaction 
density is suppressed for a hierarchical mass spectrum. 



V. RATE EQUATIONS 

Solving a system of Boltzmann equations in general re- 
quires the use of numerical codes capable of treating large 
systems of stiff differential equations for the different mo- 
mentum modes. This is a difficult task if one wants to 
study a wide range of model parameters. A commonly 
employed simplification is to approximate the Boltzmann 
equations by the corresponding system of 'rate equations' 
for the abundances Ya- In |35| it was shown that the two 
approaches, Boltzmann or the rate equations, give ap- 
proximately equal results for the final asymmetry, up to 
^10% correction. 

Starting from a quantum Boltzmann equation of the 
type (|53p we derive here the rate equation for the lepton 
asymmetry which includes the (usually neglected) quan- 
tum statistical factors. In our analysis we are closely 
following [22] . Contribution of various processes to the 
generation of the lepton asymmetry can be represented 
in the form: 



2?m/ 



"'^^Niab...jk 



E 

[-'Niab...^jk... ^Niab...<r^ij. 



r 



.gPaPb...;PjPfc...„ 
Nia.b...<-*Jk... ^Niab...<r^ij. 



(59) 



compare with (|53p . where the sum runs over each possible 
particle state with £ € {j} or i G {a}. We assume that 
the SM particles arc maintained in kinetic equilibrium 
by the fast gauge interactions. This means that their 
distribution function takes the form; 



fait,Ea) 



Tl) 



(60) 



with a time- (or temperature-) dependent chemical po- 
tential ^a = IJ^a{t)- Here the upper (lower) sign corre- 
sponds to bosons (fermions). It is also useful to define 
the equilibrium distribution function. 



re 
J a 



eq 



.Ea/T 



Tl 



(61) 



The fast SM interactions relate chemical potentials of 
the leptons, quarks and the Higgs, such that only one 
of them is independent. We can therefore express the 
chemical potential of t he q uarks as a function of the lep- 
ton chemical potential |36l - l38l | , 



1 



fJ-t 



21 



Hi = cuHe, , hq 



-Ui = cgefj-e . 



Chemical potentials of the antiparticles: /ig = —fJ-a- The 
lepton chemical potential is related to the abundance by: 



T 



ce 



2Y, 



eq 



(63) 



where ci depends on the thermal mass of the lepton. For 
mi/T « 0.2 it can be very well approximated by the zero 
mass limit, ci « 9^(3)/7r^ « 1.1. 



Using the identity 1 ± /a = e 



(£» 



.)/T 



fa and energy 



conservation we can rewrite the combinations of distri- 
bution functions appearing in (j59p as: 



-pqPaPb---;PjPk- 
•'Niab...-i^jk... 



(64) 



(2^rM9+EaPa-E,ft 






[fm a - ai^ - /A^.){e^^ '''^'^-^ '^'^ - 1}] , 



where we have suppressed the momentum argument in 
the distribution functions for notational convenience. We 
can then expand (|64[1 in the small chemical potential ^a- 
Assuming the Majorana neutrino to be close to equilib- 
rium, f^i — f'j^. ^ C(Ma), we see that the term in the 
square bracket in (|64p is already of the first order in the 
chemical potential. We can therefore replace the distri- 
bution functions in the second line of (|64p by the equi- 
librium ones, 



n f^^n fl± f"^) 



1 



req 

•f 1 



Ni 



n/rn(i±/a')' (^s) 



and expand the exponential to first order in the chemical 
potential. Since we assume small departure from equi- 
librium the Majorana distribution function that multi- 
plies the chemical potential should also be replaced by 
the equilibrium one. The corresponding equation for the 
antiparticles is obtained from the above equation by re- 
placing /ia -^ —^ia- The last step is to assume that 
the Majorana distribution function is proportional to its 
equilibrium value. 



In. 



YnA^) 



req 



Ni ■ 



(66) 



Putting everything together we get the conventional form 
of the rate equation, 



sUdYL 

z dz 



El Niab... Niab...\ I XEj. 
K^jk... "fjk... I I yeq 



-1 



Yr 



I I\iab...\ ^ 1^ 

\'jk... /(^e'^ab...^jk... riY^q 



TYt 



(67) 



(,6^j where we have defined the production and washout reac- 
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tion densities 

/ Niab... Niab... 



^jk.. 



X eN,ab...-^jk... y^N,ab...^jk... '^ ^Niab...^jk... J Jn, 

X (n fa' 11(1 ± /;«) + n /f n(i ± /a"^)) ^ (ssa) 



(7. 



Afia&...\ I 
jfc... /|w 



dIi^X.Z'-'^-f^[^ + ^aPa-E,P^ 



X 



(^"Ariab...<^jfc... + "JVia6...4^iA;...J /AfiV-*- JNi) 

(n /a' 11(1 ± /f ) + n /r n(i ± /a^) 



J 



and the numerical factor, 

_ Ej A*J - Ea A^a 



Cab...-<r^jk.. 



fit 



(68b) 



(69) 



Equation (|67p must be supplemented by an equation for 
the heavy neutrino abundance, 



sndY, 



Ni 



dz 



El Niab 



N, 



- 1 



(70) 



with the reaction density given by an expression similar 
to (|68bll : 



(7 



Niab.. 
jk... 



y^Niab...<-^jk... ~^ ^Niab...^]k...J JNi 

(n fa' 11(1 ± /;') + n z;'^ n(i ± fa 



(71) 

a J J a 

Note that these expressions are valid for two-body scat- 
tering processes with Majorana neutrino in the initial or 
final state as well as for Majorana (inverse) decay pro- 
cesses. 

If the quantum-statistical corrections are neglected, i.e. 
if the 1 ± / terms are replaced by unity and the equi- 
librium fermionic and bosonic distributions are approx- 
imated by the Maxwell-Boltzmann one, then (j68bp and 
([7T|) are equal. For the case of a 2 o 2 scattering process 
they read: 



\7jfc' ) — / ""-'-Niajk 



dnT'i,'"'{2^f5{q + pa-v,-Pk) 



(72) 



^ {^Nia^jk + ^Nia^jkjJNifa ■ 

Part of the integrations in ((7^ can be performed analyt 
ically and we obtain: 



(7f"^ 



>jk 



T 
64^ 



OC 



(73) 



Here Sr. 



{Mi + ma)"^ (assuming Mi + rua > mj+mk) 



and a-{s) is the so-called reduced cross-section: 



-Nil 



27r 
1 f dipa 



in J 2tt 
t 



dt 



{^Nia^jk + ^Nia*^jk) ' \'^) 



where s and t are the usual Mandelstam variables. The 
integration limits are given by J22l | : 

t± = Mf + m] 

- ^ [(1 + M}ls - m2/,s)(l + m)ls - mils) 

T A^ (1, Mf/s, mlls)\^ (1, m^/s, m^/s)] , (75) 

where A(a;, y, z) = a;^ + y^ + z^ — 2xy — 2xz — 2yz is the 
usual kinematical function. If effective thermal masses of 
the SM particles are neglected then the integration limits 
simplify to i+ = and t~ — —{s — Mf). Integrating 
(IHl) and p2p over t and neglecting the effective masses 
of the initial and final lepton and quar ks we obtain the 
standard expressions (see, e.g. [3^, |40|) for the reduced 
cross-sections of the Higgs-mediated scattering processes: 



-NiQ 



<y^ 



9w9s 



X — 2ai + 2a^ 



{h^h)u |A|- 



(76a) 



2a^ 



-Nil 



X — ai + t 
{h)h)u\\[' 



9^9s^u-tu\ 1x12(2^-0*) 



X — Qi 
1 



In 



a; — ai -t- 1 
0^4, 



{x- a^y 



(76b) 



where we have replaced s by a; = s/Ml and introduced 
dimensionless quantities at = Mf/Mf and a^ = rn^/Mi. 
Combined with (|73)) . expressions ((76)l give the conven- 
tional reaction densities of the Higgs-mediated scattering 
processes. 

Since in the conventional approach the CF-violating 
parameter is calculated in vacuum it is momentum- 
independent and therefore can be taken out of the in- 
tegral. The CP-violating reaction densities are thus pro- 
portional to the washout ones: 



,w>< 



,Nii 



-jk Ijk 



(nT) ■ 



(77) 



where we have again assumed a strongly hierarchical 
mass spectrum of the heavy neutrinos. When the 
medium corrections are taken into account the CP- 
violating parameter depends on the momenta of the ini- 
tial and final states and this simple relation is violated. 



VI. NUMERICAL RESULTS 

To illustrate the effect of the quantum-statistical cor- 
rections and effective thermal masses we present in this 
section ratios of the reaction densities to the conventional 
ones assuming a strongly hierarchical mass spectrum of 
the Majorana neutrinos. 
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/ Nia\\ 1 1 Nia\ 



,NiaNia'..,NiaNia^ 

y^jk Ijk /ly^jk 'jk /MB,m=0 





FIG. 8: Ratios of the scattering reaction densities obtained 
taking into account the thermal masses (dashed lines) and the 
thermal masses plus quantum-statistical effects (solid lines) 
to the conventional ones. The thick solid lines correspond to 
(|68b|l whereas the thin ones to (|71|) . 



Let us first consider the scattering processes. Ratios 
of the improved reaction densities to the conventional 
ones are presented in Fig. [5] Note that the Majorana 
(as well as the quark) Yukawa couplings cancel out in 
these ratios and for this reason we do not specify them 
here. The dashed lines show the ratio of the reaction den- 
sity computed using ([751) - (|75|) . i.e. taking into account 
the effective thermal masses but neglecting the quantum- 
statistical corrections, to the conventional ones. For the 
Nil o Qt process (dashed red line) the effective masses 
decrease the available phase space and lead to a suppres- 
sion of the reaction density in the whole range of temper- 
atures. Note that the ratio does not approach unity at 
low temperatures. Qualitatively this behavior can be un- 
derstood from ((72|) . Let us assume for a moment that the 
scattering amplitude is momentum-independent. The re- 
action density at low temperatures can then be estimated 
by evaluating the distribution functions at the average 
momenta (pi) and (pa) ~ 3T. In the ratio of the reaction 
densities the Majorana distribution function cancels out 
and 



(7y)A//a,m#o _ cxp(-£^a/r) 

(7r)MB,m=0 exp(-(pa)/T) 



exp{-ml/2{pa)T) 



A more accurate estimate for the ratio of {'~fY)MB,m^o 
and (7y )MB,m=o is ^ exp(— m^/T^). Since to a good 
approximation nia oc T we conclude that this ratio is a 
constant smaller than unity. In other words, despite the 
fact that at low temperatures the quark masses become 
small compared to the Majorana mass this ratio is not 



FIG. 9: Ratio of the CP-violating reaction densities to the 
ones computed using Boltzmann statistics and neglecting the 
thermal masses of the initial and final states. 



expected to approach unity as the temperature decreases. 
Note also that (in a very good agreement with the numer- 
ical cross-check) this ratio does not depend on the masses 
of the final states. Of course, the momentum dependence 
of the scattering amplitude somewhat changes the low- 
temperature behavior of the reaction density. Interest- 
ingly enough, for the NiQ ^ it process the inclusion of 
the thermal masses actually enhances the reaction den- 
sity at high temperatures (dashed blue line). This occurs 
because the induced increase of the amplitude turns out 
to be larger than the phase-space suppression. At low 
temperatures the effective masses become negligible in 
the scattering amplitude but still play an important role 
in the kinematics. As a result, the ratio becomes smaller 
than unity and continues to decrease as the temperature 
decreases. Let us note that for a (moderately) strong 
washout regime most of the asymmetry is typically pro- 
duced at z < 10 and the low-temperature behavior of 
the reaction densities does not affect the generation of 
the asymmetry. Since all particles in the initial and final 
states are fermions the quantum-statistical effects further 
suppress the reaction densities (solid blue and red lines) 
and render the ratio of the improved and conventional 
reaction densities smaller than unity for both NiQ o £t 
and Nil o Qt in the whole range of temperatures. Ratios 
of the improved CP-violating scattering reaction densi- 
ties to the conventional ones are presented in Fig. [21 For 
both scattering processes the improved CP-violating re- 
action densities are enhanced at high temperatures. This 
is explained by the enhancement of the CP- violation in 
the Majorana decay observed in [l^ [221 ■ ^^ the inter- 
mediate temperatures the relative enhancement of the 
CP-violating parameters gets smaller and is overcompen- 
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FIG. 10: Ratio of the Ni -f-)- IQt decay reaction density ob- 
tained taking into account effective thermal masses and quan- 
tum-statistical effects to the conventional one computed tak- 
ing into account only the effective thermal masses of the final 
and intermediate states. The thick solid line corresponds to 
whereas the thin one to H71[l . 



sated by the effective mass and Fermi-statistics induced 
suppression of the washout reaction densities that we 
have observed in Fig.[51 The low-temperature behavior is 
somewhat different for the two scattering processes. For 
the Ni(, o Qt process the effective mass and quantum- 
statistical effects get smaller in both the (unintegrated) 
CP-violating parameter and the washout reaction den- 
sity, such that the ratio of the CP- violating reaction den- 
sity to the conventional one slowly approaches a constant 
value. On the other hand, for the NiQ o tt process the 
suppression of the washout reaction density induced by 
the effective masses of the initial and final states that 
we observed in Fig.[S] also leads to a suppression of the 
CF-violating reaction density that gets stronger at low 
temperatures. 

Next we consider the three-body decay. Neglecting the 
quantum-statistical effects and using vacuum approxima- 
tion for the Ni o iQt decay amplitude in (|68p and ((7T|) 
we recover the conventional expression for the decay re- 
action density: 



\ hQtl 



i^TMf^N,- 



,,Q,K,{M,/T). (78) 



Note that it is important to retain the effective thermal 
masses of the quarks in the calculation of ^ Ni^iot- The 
four-momentum of the intermediate Higgs, see Fig. [51 
varies in the range {mq + mtY' ^ k^ ^ [Mi — miY . 
The relation m^ < mq + mt, which is fulfilled in the 
SM, ensures that the intermediate Higgs remains off-shell 
and prevents a singularity in the Higgs propagator. The 



FIG. 11; Ratio of the CP- violating reaction density of the 
Ni -H- £Qt process obtained taking into account effective ther- 
mal masses and quantum-statistical effects to the ones com- 
puted taking into account only the effective thermal masses. 



ratio of the reaction density computed taking into ac- 
count the quantum-statistical effects and effective masses 
to the one computed taking into account only the ef- 
fective masses is presented in Fig.[Tni Note that since 
ruQ pa mt ~ 0.4 T and mi « 0.2 T this three-body decay 
is kinematically allowed only at T < Mi. As one would 
expect, at high temperatures the fermionic nature of the 
initial and final states leads to a suppression as compared 
to the Boltzmann approximation. At low temperatures 
the quantum-statistical effects play no role and the ra- 
tio slowly approaches unity. Ratio of the CF-violating 
reaction density for the Ni o £Qt process is presented 
in Fig. 111! At high and intermediate temperatures the 
medium enhancement of the (unintegrated) CP- violating 
parameter is overcompensated by the suppression of the 
washout decay reaction density that we have observed 
in Fig.[Tni At low temperatures the effective mass and 
quantum-statistical effects get smaller in both the (unin- 
tegrated) CP-violating parameter and the washout reac- 
tion density, such that the CP-violating reaction density 
slowly approaches the conventional one. 

To conclude this section we present the ratio of the 
three-body decay and 2 o 2 scattering processes to the 
reaction density of N <-;• £(/) process, see Fig.[T21 As can 
be inferred from this plot, the three-body decay is sub- 
dominant in the whole range of temperatures and can be 
safely neglected. The inclusion of the effective masses has 
a very similar effect on the two-body decay and scattering 
reaction densities such that the ratios of the two almost 
do not change as compared to the one computed in the 
massless approximation. The inclusion of the quantum- 
statistical corrections has a stronger effect on the scatter- 
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FIG. 12: Ratio of the washout scattering and three-body de- 
cay reaction densities to the reaction density of Ni -o- Itj) pro- 
cess. The dashed hnes denote the ratios of the conventional 
reaction densities, the thin solid lines the ratios computed 
taking into account only the effective masses in all the reac- 
tion densities, and finally the thick solid lines the ratios com- 
puted taking into account the effective masses and quantum- 
statistical corrections in all the reaction densities. The reac- 
tion density (7f^'),„ is computed using (O and the definition 
(I68bll. see also 



ing processes sueh that the ratio of the reaction densities 
is smaller than the ratio of the conventional ones. Let us 
also note that the scattering processes are very impor- 
tant at high temperatures but become subdominant at 
low temperatures. 



tons and quarks. 

We find that the inclusion of the effective masses and 
quantum-statistical terms suppresses the washout reac- 
tion densities of the decay and scattering processes with 
respect to the conventional ones, where these effects are 
neglected, in the whole relevant range of temperatures. 
For the Nil o Qt process the ratio of the improved 
and conventional washout reaction densities slowly ap- 
proaches a constant value close to unity at low tempera- 
tures. Interestingly enough, for the NiQ <-)• it processes 
this ratio decreases even at low temperatures. Finally for 
Ni -H> £Qt process the ratio slowly approaches unity at 
low temperatures. As far as the CP- violating reaction 
densities are concerned, we find that for the scattering 
processes the ratio of the improved and the conventional 
ones is greater than unity at high temperatures but is 
smaller than unity at intermediate and low temperatures 
because of the thermal masses and quantum-statistical 
effects. For the three-body decay this ratio is smaller 
than unity in the whole relevant range of temperatures. 

We expect that the effects studied here can induce a 
0(10%) correction to the total generated asymmetry. For 
a detailed phenomenological analysis it is necessary to 
include further phenomena such as flavour effects and 
process with gauge bosons in the initial and final states. 
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VII. SUMMARY 

In this work we have studied AL = 1 decay and scat- 
tering processes mediated by the Higgs with quarks in 
the initial and final states using the formalism of non- 
equilibrium quantum field theory. 

Starting from the Kadanoff-Baym equations for the 
lepton propagator we have derived the corresponding qu- 
antum-corrected Boltzmann and rate equations for the 
total lepton asymmetry. As compared to the canonical 
ones the latter are free of the notorious double-counting 
problem and ensure that the asymmetry automatically 
vanishes in thermal equilibrium. To compute the colli- 
sion term we have taken into account one- and two-loop 
contributions to the lepton self-energy and used the ex- 
tended quasiparticle approximation for the Higgs propa- 
gator. The impact of the SM gauge interactions on the 
collision term has been approximately taken into account 
in the form of effective thermal masses of the Higgs, lep- 



Appendix A: Higgs self-energy 

The top quark contribution to the Higgs self-energy is 
derived from the 2PI effective action. At one-loop the 
contribution of the top quark is given by: 

iT2 =5s|^P / d'^ud*vtY[SQ,^^{v,u)PRSt{u,v)] 

>^4Acdiv^u)e^a: (Al) 

where the factor gs ~ 3 comes from the summation over 
color indices and e = ia2- In a SU(2)l symmetric state 
the Higgs and lepton propagators are proportional to the 
identity in the SU{2)l space, and so is the Higgs self- 
energy. 



^abi^^y) =^ix,y)Sab 



Sir. 



\aiy^x) 



gs\\\hT [Sq (y, x)PRSt (x, y)PL] 5ab ■ (A2) 
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Its Wightman components are given by, 

n^ix,y) = gs\X\hT[SQ^iy,x)PRSt;,ix,y)PL] . (A3) 

Finally, performing a Wigner transform of the above 
equation, we find. 



X ti-[SQ^it,pQ)PRSt^it,pt)PL] . (A4) 

Appendix B: Reaction density of 1 — !> 3 decay 

For Ni — > £Qt decay the general expression (j7ip takes 
the form: 



{<6t)= <^Zo:(^-rs(<i-p-pQ-p^) 



NUQt 



X ^M.^iA, X A B^ 4(100 +Pt) X ^^ 



(Bl) 



X /^. [(1 - /r)(i - /|')(i - r ) + tPfQ'ft'] , 

where we have used the explicit form of the decay ampli- 
tude p3)) . To reduce it to a form suitable for the numer- 
ical analysis we insert an identity: 



ds d kS{pQ + pt ~ k)6+{k — s) 



(B2) 



where k is the four-momentum of the intermediate Higgs. 
Approximating furthermore A^_|_^ by A^ we can rewrite 
the reaction density in the form: 



\ hot/ 






ds 

N,JN, J ^ 



A|(.) 



(B3) 



X [(1 - /r)(i - /o')(i - in + frfo'in ■ 

Note that in the regime m^ < tuq + rrit realized in the 
considered case the Higgs is always off-shell and its width 
can be neglected in A^. For the second line in (jBSp we 
can use 12211: 



JdUf^{27rfS{q^k~p) 



-> 



8^|q| 



dE, 



dip 

''' I 2^ 




The integration limits are given by 

E^ = ^[Eq{l + xe-x^)± |q|A5(l,a;^,x^)] 



(B4) 



(B5) 



where xg = m'j/M^, x^ = s/Mi and X{x,y, z) = x^+y'^ + 
z^ — 2xy — 2xz — 2yz is the usual kinematical function. 



For the third line we can use a similar expression with 
xq = riiql s and Xt = rn^/s. 

Expressed in terms of the integration variables the am- 
plitudes take the form: 



^N,^i<p = gw{h'^h)u{Mf + mj- s), 



(B6a) 
(B6b) 



Since they do not depend on the angles between the 
quarks and leptons the integration over ip is trivial and 
the reaction density takes the form: 



{iZ)-jd^kaJY.^i{^) 



(B7) 



""t dEr, 



■ ^Ni^l4> 



Kq dE, 



PQ ^ 



E-^ 8^ikr^->Q* 



IE- 8^|q| 

X [(1 - /;«)(i - n^){i - fn + frro'ft'] 



The three-momentum of the intermediate Higgs is given 
by |k| = {El - s)^ and Ek = Eg - Ep. Note that if 
we neglect the quantum-statistical factors in (jB7p the 
reaction density takes the standard form. 



Appendix C: Majorana spectral self-energy 

We compute here the Majorana spectral self-energy. In 
a CP-symmetric medium it reads J22i |: 



p 



^[{h^h),,PL + {h^h)*^Pn]Ls, (CI) 



where we have defined the loop function Ls{q), 

Ls{q) =167r jdn% {2n)H{q-p-k)f 

X [A^(fc)S,(;,)(p) + A^(,)(fc)S^(p)] . (C2) 

Using the eQP for the Higgs, see (|45p . one can split the 
function Lc into a decay part, identical to the one com- 
puted in j23 |. 

Lf (g) = 16^1 dnliFlSl^,^f , (C3) 

where we have assumed q^ > Q, and defined 



r, 



aPb---\PiPj--- _ 

{a)b...<r^ij... — 



(27r)' S{pa +Pb+ ...-Pi-Pj - ■■■) 
X 



[/r...(i±/f)(i±/f)... 
+ /f/f ...(i±/r)---] , (C4) 

and a scattering part, 

LUq) = l67rl dUl,p^p^{27T)^S{q+pQ-p-pt) 

X Sp{p)SQp{pQ)Stpipt)A'^R+A{Pt - PQ)^^i^Qp 

X [C(i - /n(i - /f) + /r/r(i - C)] • (cs) 
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Performing the frequency integration as explained above, 
see (|55|) . we can rewrite (jCSp as a sum of four terms, 
corresponding to the three scattering and one three-body 
decay process as weh as their CP-conjugates. Assuming 



L%{q)^lQ.jdI%r 



+ -f-^N^yi^Qt^R+AyPt +PQ)-4,-,Qt 



(C6) 



In the regime m^ < mq + mt the intermediate Higgs 
cannot be on-shelL Therefore one can neglect the Higgs 
width in A^_^^ and approximate it by Ay. 

As can be inferred from the definition (jC4p for the scat- 
tering terms J- vanishes in vacuum, whereas for the decay 
term it does not. Due to Lorentz covariance in vacuum 
both L'g and L| must be proportional to the four- vector 
q. Using (jB4p and (jBSp we find that the coefficient of 
proportionality is equal to unity for the decay contribu- 
tion, i.e. L'l, ~ q, if thermal masses of the Higgs and 
leptons are neglected. Using (jB7p we find that for the 
scattering contribution the coefficient of proportionality 
reads: 



above, if thermal masses of the Higgs and leptons are ne- 
glected then Lg = g in vacuum and we recover e = e™'^. 
Once the contribution of the Higgs-mediated processes is 
taken into account e"""^ ^ e™'^. To estimate the size of the 
corrections induced by (jC6p we plot the ratio of thermally 
averaged Cf-violating parameter, (e) = {^'^n) / (In) ^ ^^ 
e™'^, see Fig. [131 Note that we have neglected thermal 
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g,\Xl /•(^•^.-™^)^ ds 

' (mq+nit) 



-,« 2 / -T72^''{'^^XQ,Xt)\^{l,Xl>,X^) 



(s- 



mt 



ml){Mf + ml ^ s) 



[s-mlY 



(C7) 



Note that since we have omitted the Higgs decay width, 
this expression is convergent only if m^ < mq + mt . The 
vacuum result (jC7p provides also a very good approxima- 
tion for nonzero temperatures provided that M/T ^ 1. 
The thermal masses of the quarks then ensure that the 
Higgs remains off-shell and therefore that (|C7|) is finite. 
It is important to note that due to the temperature de- 
pendence of the effective masses the coefficient (|C7p is 
temperature-dependent as well. A numerical analysis 
shows that it grows as the temperature decreases. 

Using Ls we can calculate the in-medium CP- violating 
parameter in Ni o 1(f) process. For a hierarchical mass 
spectrum |22l |: 



, pLs 
qp 



(C8) 



where e™'^ denotes the vacuum CF-violating parameter 
calculated neglecting contributions of the Higgs-mediated 
processes, i.e. using only Lg. As has been mentioned 



FIG. 13: Ratio of the thermally averaged CP-violating pa- 
rameter to the one calculated in vacuum neglecting the con- 
tribution of the Higgs-mediated processes. The blue line cor- 
responds to (|C3|I . whereas the red lines to the sum of (|C3p 
and (|C6|I . The dashed red line is obtained by omitting the 
contribution of the three-body decay in (|C6p . 



masses of the final-state Higgs and lepton in the numer- 
ics. The blue line corresponds to the CP-violating pa- 
rameter computed using (jCSp . In agreement with the 
above discussion the ratio reaches unity at low temper- 
atures. The dashed red line corresponds to the CP- 
violating parameter computed using the sum of (|C3[) 
and the scattering (lines two to four) contributions to 
(jC6|) . As expected, at high temperatures we observe an 
enhancement of the ratio, whereas at low temperatures 
it reaches unity. The solid red line is obtained by con- 
sidering the sum of (jCSp and all of the terms in (jC6p . 
Since the three-body process is kinematically suppressed 
at high temperatures, the dashed an solid lines overlap 
for z < 1. At lower temperatures the quantum-statistical 
effects are small. However, in agreement with the dis- 
cussion below (jC7[) . the effective thermal masses of the 
Higgs and quarks lead to a slow rise of the ratio at low 
temperatures. 
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